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ABSTRACT 



Aims. The spectacular variability of microquasars has led to a long string of efforts in order to classify their observed 
behaviors in a few states. The progress made in the understanding of the Quasi-Periodic Oscillations observed in these 
objects now makes it possible to develop a new way to find order in their behavior, based on the theorized physical 
processes associated with these oscillations. This will also have the interest of reuniting microquasars in a single 
classification based on the physical processes at work and therefore independent of their specificities (mass, variation 
timescale, outburst history, etc.). This classification is aimed to be a tool to further our understanding of microquasars 
behavior and not to replace phenomenological states. 

Methods. We start by considering three instabilities that can cause accretion in the disk. We compare the conditions 
for their development, and the Quasi-Periodic Oscillations they can be expected to produce, with the spectral states in 
which these Quasi-Periodic Oscillations are observed and sometimes coexist. 

Results. From the three instabilities that we proposed to explain the three states of GRS 1915+105 we actually found 
the theoretical existence of four states. We compared those four states with observations and also how those four states 
can be seen in a model-independent fashion. Those four state can be used to find an order in microquasar observations, 
based on the properties of the Quasi-Periodic Oscillations and the physics of the associated instabilities. 
Conclusions. 

Key words. X-rays: binaries, stars: individual (GRS 1915+105, XTE J1550— 564), accretion disks 



1. Introduction 

Since the first observations of microquasars, their strong 
variability has led to successive attempts at classifying 
them in universal spectral states, defined by observables 
such as the luminosity, the energy spectrum of the emis- 
sion, and also the presence of Quasi-Periodic Oscillations 
(QPOs). Historically, the classification has involved five dif- 
ferent states, based on the shape of the energy spectrum 
and the flux level in the 1-10 keV band, (see for example 
IVan der Klisl . 11994 iNowakl . I1995D . The first classification 
attempts started with the observation of Cyg X-l and GX 
339-4, and interpretations were based on the mass accre- 
tion rate of the black hole. As M increased the source was 
thought to go from quiescent to low/hard, intermediate, 
high/soft and very high state. 
These states can be described as follows: 

Low/hard state The spectrum is a power-law with a pho- 
ton index of the order of 1.5 to 1.9 and an exponential 
cutoff around 100 keV. The X-ray luminosity is esti- 
mated to be below 10% of the Eddington luminosity 
(INowakl . 119951b 

High/soft state In this state the total luminosity is 
higher; the spectrum is dominated by a soft, blackbody- 
like component with a characteristic temperature of the 
order of 1 keV. A power-law tail is also present but much 
less luminous than in the hard state, and its photon in- 
dex is close to 2.5 



Intermediate state This state is seen during transitions 
between the low/hard and high/soft states, with spec- 
tral characteristics ranging 'in-between' those of these 
states. 

Very high state In some systems a state with very high 
luminosities is observed. In this state the nonthermal 
tail and blackbody components become comparable in 
flux. The power-law component has a photon index of 
2.5 and does not show evidence of a cutoff even out to 
a few hundred keV. 

Quiescent/off state In the last state the object appears 
to be 'off', with a flux level much lower than in the other 
states and a photon index softer than in the low/hard 
state. 



iRemillard et al\ (|2002f) . and again lRemillard fc McClintocld 
(|2006l ) took another approach. They defined a set of three 
'fundamental' states (using more descriptive names, 
not based on the X-ray luminosity) and the transitions 
between them. These three states are the Low-Hard, Steep 
Power Law (Very High) and thermal (High-Soft) states, 
respectively LHS, SPL and thermal states hereafter. Other 
states are simply transitions between these three. XTE 
J1550-564 is a good example of this classification. 
Given the strong differences between the behaviors of the 
known microquasars, the differences between these classifi- 
cations can arise simply from the sources that receive more 
attention in each work. It is however interesting to note 
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already the strong link between these three classes and the 
temporal behavior of the sources: the thermal state does 
not show QPOs and the SPL has both low-frequency and 
high-frequency QPOs (hereafter LFQPO and HFQPO). 
For the low-hard state the association is not as straightfor- 
ward. This state often exhibits a LFQPO whose frequency 
has been shown to be correlated with various properties 
of the accretion disk, and therefore ultimately to the 
accretion rate itself. One must note, however, that there 
are sources where the LFQPO is not detected in the 
LHS. There are several possibilities for this absence. 
Either the QPO is truly absent, or the amplitude of the 
QPO is too weak to be detected. The latter possibility 
can have various origin again. The low luminosity of the 
source may not permit to detect it, or the QPO frequency 
might simply fall out of the frequency range. Another 
explanation is that the underlying instability mechanism 
giving rise to it is itself weak or fails to generate a strong 
observable modulation of the X-ray emission. This can 
be the case if, e.g., the physical properties of the disk 
and its corona, or of its inclination (one must remember 
for instance that a simple rotating hot spot would not 
cause a QPO, even in a strongly inclined source) make the 
detection difficult, or the QPO falls out of the frequency 
range analyzed, or the QPO is truly absent. In a weak 
source such as XTE J1752-322, the LFQPO seems absent 
from the LHS, but appears weakly when observations are 
summed ( Munoz- Daria et al. I. 201d[). In a stron ger source 
such as XTE J1550-564 (jRodriguez et aai2003D . LFQPOs 
appear early in the initial Low-Hard state and immediately 
after the transition back to the Low-Hard state. Given the 
very frequent observation of LFQPOs in the LHS in other 
sources (H1743-322, e.g. Prat et al. 2009; XTE J1118+480, 
Revnivtsev et al. 2000; XTE J1908+094, Gogus et al. 
2004; 4U 1630-474 Dieters et al 2000) we consider it as a 
property, which may or not grow to be observable, of that 
state. 

The original three state description has been refined 
further using from the start X-ray timing properties, 
especially the type of Low-frequency Quasi-Periodic 
Oscillation, and observations at other wavelengths such 
as the presence of a jet seen in ra dio, t o separate the 
different states. iHoman fe Belloiiil (|2005h introduced a 
4-state classification in the continuity of the previous 
ones, but with a stronger link to the QPOs observed, 
and in particular to the identification of three different 
flavor s (A, B and C) of the LFQPO (jRemillard et all 
120021 ). The High-Soft state is in the direct continuity of 
the previous classification while the Low-Hard state is 
defined similarly with the exception of the observation 
where a type C LFQPO is in the observation. It mainly 
correspond the early part of outbursts. The two other 
states of their classification are the Hard Intermediate 
State, in which a C-type LFQPO is detected, while the 
Soft Intermediat e state is characterized b y either a B- or 
A-type LFQPO (|Homan fc Bellonil . I2005T) . The differences 
between those LFQPOs include the evolution of both 
frequency and amplitude, the harmonic content and the 
noise components such as the Band Limited Noise. This 
classification is coherent with an evolution from Low/Hard 
to Hard Intermediate, Soft Intermediate, High Soft, Hard 
Intermediate and back to Quiescent/Low-Hard State. 
However, even this more refined classification does not 



encompass all the observed behavior; for example the case 
of HFQPOs without LFQPOs (as seen in GRS 1915+105 
in the B state) does not fit in any class. 



A different classification was defined by iBelloni et al\ 
(|2000l ) for the special case of the most spectacular and 
active microquasar, GRS 1915+105. They defined 3 funda- 
mental spectral states, labeled A, B and C (not related to 
the A, B and C LFQPO types), which combine into 12 re- 
producible classes of variability (a, /3, ...). The three states, 
defined in a color-color diagram, are again characterized by 
different contributions from the disk (thermal component 
of the X-ray) and its corona (non-thermal component at 
higher energy). Each of the 12 variability classes is either 
composed of only one state A (class </>) , C (x) or repetitive 
cycles of various length of BAB (6, 7), CAC (9) or CAB 
(a, v, p, k, A, P, fi). This classification has proven useful, but 
is limited to that particular object. The main difference 
between GRS 1915+105 and the other microquasars, 
besides the fact that it has remained active since it was 
first observed, is actually the occurrence of these cycles on 
short (up to a few tens of minutes) timescales, and their 
repeatability over the years. 

There is no one-to-one connection between this spectral 
classification and the previous ones, partly because GRS 
1915+105 tends to sta y at a much high e r flux level than 
other objects (see e.g. iFender fc Bellonil [20041 ). However, 
considering the short-term temporal behavior of the three 
basic spectral states, a one-to- one associatio n does seem 
to exist: indeed, the A state of IBelloni et all (|2000h does 
not show any quasi-periodic oscillation, while C always 
exhibits the Low-Frequency one. The B state seems to 
always be present when HFQPO are detected alone (i.e. 
without a LFQPO) though no extensive study has been 
published yet (see for example Morgan et al., 1997). 

These classifications are based on observations and 
their differences point to the differences between the 
sources used as reference. The SPL, which might actually 
be defined by the joint observation of LF and HFQPOs, 
was often seen in XTE J1550-564 but was not, until 
recently, observed in GRS 1915 +105 and th e refore was 
absent from the classification of IBelloni et al\ (|2000j). As 
we have alread y mentioned, t he sam e can be said about 
the B state of IBelloni et al\ (|2000f ). during which GRS 
1915+105 seems to exhibit a HFQPO alone. This has 
not yet been observed in XTE Jl 550-564 and is therefore 
abs ent from the classification of Remillard et ai\ (|2002f ) 
and iRemillard fc McClintockl (|2006l) . 

Here we take another approach. We decided to look at 
the different behaviors of microquasars from the perspective 
of the disk instabilities happening in them. We are focus- 
ing here on instabilities that have clear observational sig- 
natures, namely the ones that could explain the QPOs (HF 
and LF). We are making the assumption that the mech- 
anisms (instabilities) at the origin of the QPOs become, 
when they happen, the dominant effect in the disk. Also, 
these instabilities will be impacting the disk even before 
the observable "QPO" can be detected with our present ca- 
pabilities (see for example the numerical work of Caunt & 
Tagger, 2001, for the evolution of a disk with the Accretion- 
Ejection Instability). The fact that, in all the existing clas- 
sifications, the different states could always be linked with 
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a change in the QPO-content tend to agree with that as- 
sumption. Our goal is to classify the behavior using what 
occurs in the disk. We therefore start from the different 
instabilities thought to be active in the accretion disk and 
look at the different observables they would lead to. Using 
this to create a classification allows it to go beyond the 
'not yet observed' arguments as it is not based on any set 
of observations but on more general principles. 

The first step in this dire ction was m a de w ith the 
magnetic flood scenario (MFS: iTaeeer et all . 12001 which 
was presented as a possible explanation for the 30 minute 
cycle (class /?) of GRS 1915+105. This scenario started 
from the tentative association of the LFQPO with the 
Accretion- Ejec ti on In st ability (AEhlTagger fe Pellatl . 119991: 



states labelled A, B and C, based on their position in a 
color-color diagram (see figure Fig[l}. With the hardness 
ratio defined as HRI = B/A and HR2 = C/A where A 
corresponds to the channel 0—13 (2 — 5 keV), B to 14 — 35 
(5 - 13keV) and C to 36 - 255 (13 - 60keV).' 



IVarniere etal\ . l2002t iRodriguez et all . 120021 )7 and pro 



ceeded by a comparison between the observed properties 
of the source in the various states of iBelloni et al\ (|2000l ) 
and the physical requirements for and expected conse- 
quences of the developme nt of the AEI. The second step 
l|Tagger fc Varnie"rel . l2006f ) was to propose the Rossby-Wave 
Instability (RWI) as a possible mechanism for the High- 
Frequency QPO in microquasars. In these works the prop- 
erties of the disk expected to result from these instabilities 
were part of the argument for their identification with the 
QPOs. 

In the present work we continue this effort by undertaking a 
classification that starts from the physics of the different in- 
stabilities that can exist in disks. Throughout this paper 
we will use one instability for the LFQPOs, with- 
out entering into their different types which will be 
treated in more detail in a forthcoming publication 
(jVarniere et al\ . l2010h . and one instability for the 
HFQPOs. Comparing the conditions required for them to 
appear, and their expected consequences on the disk, we 
will then turn to the observations and seek the correspon- 
dence that can exist. As this is based on the instabilities at 
the origin of the QPOs, we will focus first on them and not 
seek any one-to-one association with the observation-based 
classifications. 

In section 2 we will thus briefly review the three states of 
IBelloni et all (|2000f ). and the three MHD instabilities that 
have been predicted in magnetized disks. We will then dis- 
cuss how these instabilities can be accommodated within 
the Magnetic Flood Scenario and the three spectral states 
it involves. In the following section we will then show that 
the occurrence of the three instabilities can actually give 
four instability-defined states. In section 4 we will thus 
concentrate on the fourth state: we will first show it to occur 
in MHD disk simulations, and find confirmation in obser- 
vations of several sources and not just GRS 1915+105 from 
which we started. In section 5 we will look a this tentative 
classification through the instabilities occurring in the disk 
and see how it can be seen as a model independent classifi- 
cation just based on physical processes without naming the 
instabilities. This will allow us to focus on what this type 
of physics-based classification can give us compared with 
the phenomenological one. 



2. Magnetic Flood: a three instabilities scenario 

2.1. The three states of lBelloni et~aT1 \200C\) 

IBelloni et al\ (|2000l ) classified the behavior of the mi- 
croquasar GRS 1915+105 in three fundamental spectral 



HRI 




HR2 



Fig. 1. three states of Belloni shown on the color-color di- 
agram of the class (5 of GRS 1915+105 . 



Their definitions of the three states, using color ratios 
HRI and HR2, are as follows: 



state A: This state is at low count rate, low HRI and 
HR2 ratios, in the lower left corner of the color-color 
diagram. In this state the disk flux represents a 
substantial contribution to the total flux. There is no 
quasi-periodic variability and the X-ray emission is 
mostly from the thermal disk emission. 
In comparison with other classifications, this is closest 
to the High/Soft or Thermal state since there is a 
dominant disk and no quasi-periodic variability. 

state C: State C has a low count rate, low HRI and 
variable HR2 depending on the length of the event. 
The power-law flux dominates with a little or no visible 
disk flux contribution. 

This state is known to exhibit a LF-QPO and Band- 
Limited Noise (BLN). 

In comparison with other classifications, this is closest 
to the Low/Hard state 

state B: This state is located above state A in the CC 
with a high HRI. This state also has a high count rate. 
Substantial red-noise variability is visible on time scales 
longer than Is. The disk is visible and hotter than in 
state A, and spectral fits indicate that it is close to 
the last stable orbit. High- Frequency QPOs alone are 
observed in this state. 

There seems to exist no equivalent to this state in the 
other classifications. 
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2.2. The three instabilities at the origin of the three states where 



Three instabilities have been proposed to exist in mag- 
netized accretion disks. They differ by their physics, by 
the conditions for their existence and by their expected 
consequences: 

Magneto-Rotationnal Instability (MRI) 

The MRI is a local instability that develops in a weakly 
magnetized disk (below equipartition with gas pressure, 
i.e. a plasma j3 = 8irp/B 2 > 1), when 



r 
t—1 



dr 



< 



(1) 



which is always true for (pseudo-)Keplerian disks. 
Numerical simulations have confirmed that it results in 
small-scale turbulence, and it is believed to be responsible 
for accretion in general. However no simulation of the MRI 
has yet shown a QPO of any kind. We thus assume that 
the corresponding variability, lacking any strong feature, 
is the one best ad apted to correspond to a power-law 
temporal spectrum (|Balbus fe Hawlevl 120021) . 

Accretion- Ejection Instability (AEI) 



The AEI ([Tagger fc Pellatl . 1 19991) is a global instability 
occurring in disks threaded by a poloidal magnetic field 
close to the equipartition, namely when the plasma /3 ~ 1. 
These conditions are the ones found necessa ry in MHD 
models of jets (see e.g. ICasse fc Ferreiral . l2000l) 

The AEI belongs to the same family as the spiral insta- 
bility of self-gravitating disk galaxies, and the Papaloizou- 
Pringle Instabil i ty. In the language of diskoseismology 
([Wagoner! . 119991 : iKatol . l200lh these could be described as 
unstable p-modes, while the RWI discussed below would 
correspond to the g-mode. The AEI is characterized by a 
spiral structure that develops in the inner region of the 
disk. At the corotation radius, where the gas and the spiral 
wave rotate at the same angular velocity, the spiral excites 
a Rossby vortex to which it transfers the energy and angu- 
lar momentum it extracts from the disk. 
Furthermore, in the presence of a low density corona, 
the Rossby vortex twists the footpoints of the magnetic 
field lines, generating an Alfven wave propagating to the 
cor ona where it might provid e the source for a wind or a 
jet (| Varniere fc Tagged. 12002] ) . This instability is stronger 
when the plasma j3 is of the order of one and requires 



7T C b >0 
or 



where 



C B = 



2f2B 2 ' 



(2) 



(3) 



f2 and n are the rotation and epicyclic frequencies (in a 
Keplerian disk f2 = k), £ is the surface density and B 
is the equilibrium magnetic field. This criterion is obeyed 
in disks with 'reasonable' density and magnetic field pro- 
files, whereas for the Papaloizou-Pringle Instability in an 
unmagnetized disk the criterion becomes 



d 



2ft£ 



(5) 



which would require a very steep density gradient to be 
destabilizing. 

We have presented the AEI as a good candidate for the 
low-f r equency quasi- p eriodi c oscillation (iRodr i guez et all 



120021 : IVarniere et~al\. 120021: iMikles et all 120091) and dis 



cussed ( Varniere fc Taggen . 20021 ) how the correlation 



(4) 



observed between this QPO and coronal activity could be 
explained in relation with the emission of Alfven wa ves. 
In a forthcoming publication (IVarniere et all 120101 ) we 
will also show how the three types of LFQPOs can be 
explained by the AEI within the framework presented here. 

Rossby Wave Instability (RWI) 

The RWI can develop in unmagnetized as well as mag- 
netized disks, requiring respectively an extremum of Cp 
or Cb- Such an extremum can be due to an extrem um of 
density (|Lovelace et all 119991: iTagger fc Melial . l2006f) or to 
relativistic effects near the Last Stable Orbit (LSO) in th e 
accretion disk of a microquasar (jTagger fc Varnierd . [20061) : 
indeed the LSO is defined by a vanishing epicyclic frequency 
k, whereas further out in the disk n is close to CI. Thus 
K 2 /2ft, and consequently Cp and Cb, have extrema near 
the LSO. This extremum of n occurs at r ~ 1.3 r^so, so 
that the RWI can occur whenever the inner radius of the 
disk is within that radius. 

As with the AEI, the RWI excites Alfven waves and can be 
expected to energize the corona if the disk is threaded by 
a poloidal magnetic field, though the field is not necessary 
for the instability. We have proposed the RWI as a possi- 
ble explanation for t he ob served HFQPOs in microquasars 
(jTagger fc Varnierel [2006]) . 

2.3. Magnetic Flood Scenario 

The Magnetic Flood Scenario (Tagger et al. 2004) was first 
i ntroduced to expla in the /? class of the classification of 
(jBelloni et all l2000t ) (also known as the 30-minutes cycle) . 
The MFS explains the repetitive X-ray behavior as a limit 
cycle determined by the advection of poloidal magnetic flux 
to the inner disk and its destruction via magnetic recon- 
nection (which can lead to relativistic ejection) with the 
magnetic flux trapped close to the source. 
The MFS starts with the identification of the Low 
Frequency Quasi- Periodic Oscillation with the AEI. It then 
proceeds by assuming that the onset of the state C of the 
Belloni et al. 2000 classification (possibly identified with the 
low-hard state in other classifications or the hard intermedi- 
ate state in the classification by Homan & Belloni, 2005.), 
which also corresponds to the onset of the QPO, is trig- 
gered when the disk magnetization becomes large enough 
(of the order of equipartition with the gas pressure, so that 
we might call this a 'fully magnetized' disk) for the AEI to 
become unstable. 

We also noted that state A was the one closest to what is 
expected from an a-disk dominated by viscous-like trans- 
port of angular momentum due to small-scale turbulence. 
We thus proposed to associate state A with the presence 
of the Magneto- Rotational Instability (MRI). Finally, we 
proposed recently the Rossby- Wave Instability (RWI) as 
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a possib le explanation for the HF QPO observed in micro- 
quasars dTagger fc Varnie"rel.l2006t) . An interpretation of the 
B state of lBelloni et all (|2000l ) was thus called for, in which 
this state would be dominated by the RWI when the disk 
inner radius is close to the last stable orbit. 
As a co nsequence we assoc iated the three fundamental 
states of lBelloni et al\ (|2000t ) with three distinct instabili- 
ties, based on their properties and variabilities. 

3. Three Instabilities but four states 

Let us now discuss the instability criteria, for the three 
instabilities of section !2.2l in terms of two parameters of the 
inner disk region: the location of its inner radius, measured 
by a parameter £ = ri nt /rpsOi an d the magnetization of 
its inner region, measured by ft = 8irp/B 2 : 

- the MRI requires a weakly magnetized disk, i.e. j3 > 1. 

It does not depend on £, and it is not believed to cause 
QPOs. 

- the RWI requires an inner edge close to the LSO, i.e. £ < 

t^ext where £ ex t is the position of the extremum of either 
Cp or Cp. It is stronger if (3 is of the order of unity, 
but can exist even in an unmagnetized disk. It has been 
proposed to explain the HFQPO. 

- the AEI requires a positive gradient of the quantity C p • 

which as discussed in section l2~2l is obtained with 'rea- 
sonable' assumptions on the radial profiles of £ and B. 
It also requires a magnetic field of the order of equipar- 
tition, i.e. ft ~ 1. Its existence does not depend on £, 
but its frequency is a fraction of the Keplerian rotation 
frequency at r.;„ t , and thus a function of £. It has been 
proposed to explain the LFQPO. 

These criteria are not mutually exclusive. In particular, if 
the inner disk edge is close enough to the Last Stable Orbit 
(so that the RWI exists) and the magnetic field is sufficient, 
the AEI can exist further out in the disk, with a lower 
frequency. 

Thus with the two parameters /3 and £ we can map four 
cases, defining four states of variability: 

P > 1) £,int > S,ext- The disk is weakly magnetized and is 
not at the Last Stable Orbit even if it is really close to 
it. The MRI must dominate, leading to a turbulent disk 
with no QPO. This state is similar to the High/Soft, 
Thermal or A state in previous classifications as it does 
not display any QPOs. 

P > 1) £,int < S,ext- The disk is weakly magnetized and its 
inner edge is close to the Last Stable Orbit. The RWI 
is present at the inner edge, while the MRI acts further 
out. Observationally this should appear as a warm disk 
with a HFQPO and no LFQPO. 

/? ~ 1, £i n t > S,ext- The disk is 'fully magnetized', i.e. the 
field is of the order of equipartition with the gas pres- 
sure, and the inner edge does not approach the Last 
Stable Orbit. The AEI causes a LFQPO so that the 
inner region of the disk cools down while the corona 
becomes active. 

/3 ~ 1, £i n t < S,ext- The disk is fully magnetized and the in- 
ner edge is close to the Last Stable Orbit. Both the AEI 
and the magnetized version of the RWI are present, pro- 
ducing both a LF- and HF-QPO. The frequency of the 
LFQPO varies little since it is a fraction of the rotation 
frequency at the inner disk edge, which cannot change 



much. The interaction between these two modes can be 
presumed to affect their characteristics and their effect 
on the disk (see Varniere et al, 2010 for a comparison 
of those effects with the different types of LFQPO). 

Thus from the two parameters j3 and £ we have defined four 
regions of parameter space, leading to four types of behav- 
ior. The three states that were discussed in the Magnetic 
Flood Sc enario are easily ass ociated with the states A, B 
and C of lBelloni et~ai\ (|2000h . Here, we arrived to the con- 
clusion that a fourth state is required to describe all possible 
behaviors. In this state the AEI and the RWI are active, 
and thus there are both the LF- and HF-QPO. 

It is also interesting to note that, while this separation 
in four state started from the QPOs alone, the presence of 
the Band-Limited Noise seems to correlate with only one 
state, namely the AEI/LFQPO-only state. 



4. The State {AEI + MRI}: /3 ~ 1, r int /r LSO ~ 1 

That fourth state is defined in the parameter space by the 
co-existence of two instabilities which lead to the observ- 
able that are the HFQPO and LFQPO. In order to validate 
this, we need to show that a state with these observables 
is actually observed and that these two instabilities can 
actually co-exist. 

4.1. Observations of a state with both HF and LF QPOs 

The MFS, based on observations of GRS 1915+105 , is at 
the origin of the association of the three instabilities with 
the three A,B and C states. The fact that the fourth state 
(the one with both LF and HF QP O) was recently observed 
in this source iBelloni et al\ (|2006l ) is important as it gives 
strength and confidence in our tentative classification. 

On the other hand, while this fourth prescription had 
already been observed in other sources such as XTE J 1550- 
564 (so this state is already part of the standard three/five- 
states description), the HFQPO-only state has so far never 
been seen in any of those other objects. There are a couple 
of possibilities for why that state has not been observed 
yet: 1) the HFQPO is too weak to be detected or 2) those 
sources are able to meet the criteria for the RWI only when 
they also meet the criteria for the AEI. In either case, 
it would be interesting to see what is singling out GRS 
1915+105 to either have its HFQPO alone strong enough 
to be detected or being able to meet the criteria for one of 
the HFQPO and not for the LFQPO at the same time. This 
would require a more extensive study of the observation and 
is beyond the scope of this paper. 

4.2. Numerical Evidence of the {AEI + RWI} state 

Since we have discussed the importance of observing the 
LF- and HF-QPO together, we now proceed to show that 
this can occur in numerical simulations too. This is numeri- 
cally demanding since it requires coping with different time 
scales and different instability criteria, and we will thus 
stick to a proof-of-principle simulation, choosing parame- 
ters that make it easier rather than claiming realism. 
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4.2.1. Numerical Setup 



As in iTagger fc Varnierd (120061) we use the 2 5D MHD 
code first introduced by Caunt fc Taggerl (| 2 lh to study 
the non-linear evolution of the AEI. The code, based on 
a Zeus-type scheme, is similar to those used for simula- 
tions of galactic spiral structures before full 3D simulations 
were possible; it describes an infinitely thin disk in vac- 
uum. The magnetic field can be described by a magnetic 
potential outside the disk, related to the field in the disk 
by a Poisson equation similar to that for self-gravity. The 
code uses cylindrical coordinates with a logarithmic radial 
grid. This allows high precision in the inner region of the 
disk where the instabilities develop, and a large dynamic 
range in r (typically 50 inner radii) to avoid boundary 
condition issues. The code also implements the FARGO 
scheme (Masset, 2000) which enhances execution speed by 
eliminating the Keplerian speed from the Courant condi- 
tion. Finally the code mimics the relativistic rotation curve 
(needed for the RW I ) with a pseudo-newtonian potential 
(|Paczvnskv fc Wiital . [198(1 : 



$ = _. 



GM 



with 



2GM/c 2 



This gives us kr = f2_R-y/(r — 3r s )/(r — r s ), which has 
a maximun close to the last stable orbit. 



4.2.2. Results 

Figure [2] shows a contour plot of the amplitude of m = 1 
and m = 2 perturbations in radial velocity, as a function 
of radius and time. These plots are a convenient way to 
study waves as they allow us to identify their propagation, 
i.e. oblique features showing a traveling wave while hori- 
zontal ones indicate a standing pattern. The AEI appears 
as a standing pattern within its corotation radius, emitting 
a wave traveling outward beyond it. The RWI is a stand- 
ing pattern trapped in the extremum of Cb and a traveling 
wave elsewhere. 

On Figf5] we see at first the m = 1 and m = 2 modea3 of 
the RWI close to the relativistic maximum of £g (which 
occurs in our simulation around 1.3 r^so)- Further away, 
we also see the m = 1 and in = 2 modes of the AEI. 




The conditions to obtain both instabilities are: 

- a disk near equipartition between the gas pressure and 
the magnetic field, namely /? = 8irp/B 2 ~ 1 

- with its inner edge near the last stable orbit, ^„ t sa 1. 

- an extremum of C B = k 2 / (2ft)£ / B 2 for the RWI 

- and a positive derivative for Cb for the AEI 



By using a pseudo-newtonian potential we already have 
an extremum of Cb close to the last stable orbit. We need 
to ensure that it also has a positive derivative further away 
in the disk so that the AEI can develop. 

We thus seek (k 2 /(20))E/B 2 to have a positive gradient 
away from the last stable orbit. For this we can use the 
Keplerian form of Cg = (f^ /2)E/i? 2 . Away from the last 
stable orbit Cg ~ Cb- We will impose Cg oc £ 01 and we 
will use it to fix the radial profile of B for a given density 
profile. This will ensure that we have the AEI developing 
in the disk. 

We then need to define the density E. Since we want a 
profile that is not too steep at the inner edge of the simu- 
lation, we take 



Fig. 2. Contour plot of modes m = 1 (a) and m — 2 (b) of 
the radial velocity as a function of radius and time. 

Having both instabilities present tends to make the 
spectrum messier compared to having them separately, as 
they are not ful l y rad ially separated. In a companion paper 
(jVarniere et al.\ . \201(i\ we will focus on how the instabilities 
are modified in that case and how it could explain some of 
the features of the different types of LFQPOs. We insist 
that, given the crudeness of the model, this numerical 
simulation can only be considered as a proof-of-principle. 
In particular the value of m that is most unstable, and 
thus dominant, depends on the profiles we have used, but 
also on the quality of t he model. Indeed in separate work 
on unmag netized disks (|Meheut et all l2010t ) we find that 
3D simulations can give a result differing from the 2D one. 



Finally we use the definition of Cg to set the magnetic 
geld: 5. An Instability-Based Classification 

I — — — 5.1. A model Independent look at the classification 

Although we have discussed it in view of definite Q 
els, the classification we propose is only based on the pre- 



" ~ y £Kq k Although we have discussed it in view of definite QPO mod- 



This setup will ensure that we have both the RWI and the 
AEI appearing in the inner region of the disk. 



1 We choose not the show all the modes but to focus of the 
first two for simplicity. 
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dieted presence or absence of LF and HF QPOs. As 

such it is model-independent and can be compared with 
other classifications starting from the spectral properties. 
Implicitly this means that we consider the QPOs to be the 
cause of the spectral properties, assumed to depend on the 
transport and deposition in the disk or the corona of the ac- 
cretion energy. Spectral-based classifications, on the other 
hand, simply describe the QPOs as properties of the spec- 
tral states. We find that these spectral properties agree well 
with what we can expect from our QPO models, but this 
may not be unique to these models and therefore any model 
agreeing with the spectral properties and having the pos- 
sibilities of both HF and LF QPOs would give the same 
four-states classification. 

This leads us to four states, defined by the presence or the 
absence of HF and/or LF QPOs. We can represent these 
four states as a tree diagram, shown in figure [3J 
This classification is actually close to the ones of 



all microquasars states 



pothesis of a strong causal link between timing and spectral 
changes. This classification aims at creating a framework to 
better understand what drives the source behavior and is 
not to be used instead of a spectral classification, easier to 
use when dealing with observation. Indeed, having a clas- 
sification based on the detection/non-detection of HF-LF 
QPOs has a limited interest when one is looking purely at 
classifying a state for reference not only because HFQPOs 
are hard to detect but also they require some processing of 
the data. 

On the other hand, we can turn around that difficulty 
and use this classification as a tool to probe for HFQPOs. 
Indeed, one can focus the search in the states that are 
clearly not-MRI or AEI- d omina ted. This is similar to the 
method iRemillard et all (|2006f) used to de tect HFQPO 
i n H1 743-322 following their discovery by iHoman et all 
(l2005h . 

In a second paper (jVarniere et all I2010D we will focus 
on the state where both the AEI and the RWI are present in 
the disk. In particular, we will study the impact of having 
both instabilities, more precisely the link between HFQPO 
and LFQPO, and if the effect on the LFQPO can be used 
to infer the presence of a weaker HFQPO. 




no LFQPOs 



HFQPOs 



no HFQPOs 



Fig. 3. Model independant view of our classification in a tree 
form. 



both iBelloni et all (l2000D and of (jRemillard et all 120061: 
IRemillard fc McClintockl 120061). The ma in difference with 
the classification of Belloni et all (|2000h . based on GRS 
1915+105, is that we separate their C state in two (AEI- 
dominated state and AEI-RWI-dominated state), depend- 
ing on the high frequency variability. Since this classifica- 
tion was first proposed, we have n ow learned that HFQ POs 
occur in the C state as defined by IBelloni et al\ (|2000t ) (see 
IBelloni et all l2006h . so that the AEI-MRI-dominated state 
has been observed in GRS 1915+105. 

Concerni ng the classification, mostly based on XTE J1550- 
564 , in (jRemillard et all . 120061: IRemillard fc: McClintockl . 
120061) . the main difference is the state with only HFQPOs 
which has not been observed yet and was not part of their 
classification. It will be interesting to look for reasons why 
that state is observed in the case of GRS 1915+105 and 
not for other objects. The major differences between our 
classification and the one of Homan & Belloni is that they 
do not have the HFQPO alone state and that they separate 
the Low-Hard state and the Hard intermediate state based 
on multi-wavelength properties (nevertheless, they clearly 
state that their components are clearly related). 

5.2. A method to find more HFQPOs? 

As one can see this instability-based classification is coher- 
ent with the other classifications based on a combination 
of spectral and timing behavior. This strengthens our hy- 



5.3. When the instability criteria are met or not 

As we showed in sect. 5.1 this model-independant classifica- 
tion in four states is only based on the presence or absence 
of the LF and HF QPO and not on the actual instabilities 
at their origin. It is interesting to notice that most objects 
do not seem to exhibit one of the four states, namely the 
one with the HFQPO-only. Up to now, GRS 1915+105 is 
the only source known to exhibit all of the four states with 
some regularity. As this source was always singled out it is 
not a surprise to find in it behaviors that other sources do 
not exhibit, but with this classification we have a way to 
tackle the problem. Indeed, there are only 2 reasons why 
those HFQPOs are not detected in other sources: 

1) the HFQPOs are present but are not detected because 
they are too weak. 

2) the instability criteria for the instability at the origin of 
the HFQPO alone (whichever instability is the 'right'one) 
is not fulfilled without the criteria for the LFQO being also 
fulfilled. 

In either case, it would be interesting to understand what 
is singling out GRS 1915+105 to either have its HFQPO 
alone strong enough to be detected or be able to meet the 
criteria for a HFQPO without that for the LFQPO. To 
answer that question would require an extensive study of 
the observation to decide in which case the other sources 
are and is beyond the scope of this paper. 

6. Conclusion 

In this paper we have discussed several disk instabilities 
thought to occur in microquasars. From their instability 
criteria we can divide the observed parameter space in four 
states characterized by their QPO content. This form a clas- 
sification of all microquasars observations as shown on the 
tree diagram^ We find that this new classification re-unites 
previous ones, which started from the spectral properties of 
the sources but, as a consequence, did not apply uniformly 
to all microquasars. 

Our discussion is guided by the instabilities that may be 
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at the origin of the variability. Following previous works we 
have taken the Accretion-Ejection Instability to be at the 
origin of the LFQPO, and the Rossby Wave Instability to be 
at the origin of the HFQPO, while the turbulent, thermal, 
disk state is assumed to be dominated by the Magneto- 
Rotationnal Instability. 

Considering the conditions for these instabilities to develop 
we find, and confirm by numerical simulations as well as in 
recent observations, that four states can be defined rather 
than the usual two, three or five states classifications used 
up to now, which can now be included in the one we present 
here. The four states are defined by no QPO, only the 
LFQPO, only the HFQPO, and both LF- and HF-QPO, 
thus covering in a simple manner all the possible observa- 
tions of QPOs. Furthermore this classification can be used 
as a guide in searching for occurrences of the HF-QPO, 
which is often difficult to detect. 
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